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1 Introduction 



Studies of bound-state properties in Quantum Electrodynamics (QED) are im- 
portant on their own for a variety of practical applications, but also stimulate de- 
velopment of theoretical tools useful in many other areas of physics. High-precision 
QED calculations, including evaluation of high orders in the perturbation theory, are 
necessary to match the precision of atomic physics experiments, particularly in the 
spectroscopy of simple atoms. In those cases where theory and experiment can both 
attain high accuracy, comparisons of measurements with predictions often enable 
determination of various fundamental physical constants, such as the fine structure 
constant, masses or mass ratios of the electron, muon, and proton, ratios of various 
magnetic moments, proton charge radius, etc. 

In other cases, when the relevant constants are known from other sources, one 
can very precisely test theoretical understanding of bound states. An impressive 
example, discussed at the previous RADCOR meeting |J, is the hyperfine splitting 
in the positronium ground state. There, the two-loop effects change the leading order 
prediction by about 12 MHz, or only 0.006%, but this effect is still larger than the 
experimental error by more than an order of magnitude! Rarely is computing of high- 
order loop effects more rewarding than it is in the spectroscopy of simple atoms like 
the positronium or muonium. 

Such high-precision comparisons of theory and experiment are possible in simple 
atoms because on the one hand, measurements can be made very accurately with 
the modern spectroscopic methods, and on the other there are hardly any principal 
limitations of the theory. Since electrons are much lighter than any hadrons, the 
computations are not hindered by non-perturbative QCD uncertainties. Spectra and 
lifetimes of simple atoms can, in principle, be evaluated with any accuracy required 
by current experiments within pure QED. 

Theoretical tools developed in this way, such as the computational techniques for 
the Feynman integrals or the machinery of non-relativistic effective theories [0, can 
subsequently be applied to solve problems in other areas of physics, such as hadronic 
properties or thermal field theory. 

In this talk I would like to present some results obtained in the last couple of 
years since the previous RADCOR meeting in Barcelona. Those results include the 
positronium lifetime, bound electron gyromagnetic factor, and a new approach to 
computing properties of bound states consisting of particles with widely different 
masses. 
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2 Positronium decay 

Lifetimes of the singlet and triplet positronium ground-states, p-Ps and o-Ps, can 
be measured with high precision ||. For several years there was a very significant 
discrepancy between the theoretical predictions for the o-Ps lifetime and the exper- 
imental results. More recently, studies performed at the University of Tokyo gave 
results consistent with the theory. The most accurate experimental results are sum- 
marized in Table [I]. Recently, the results obtained in gases were reexamined and new 
corrections were taken into account. The preliminary updated central value for the 
o-Ps decay rate measured in this method is approximately 7.047/yUS 0, in very good 
agreement with the vacuum measurement, but still significantly higher than the QED 
prediction. 

Table 1: Recent experimental results for the o-Ps lifetime. "Method" in the second column 
refers to the medium in which the o-Ps decays. The last column shows the value of the 
two- loop coefficient B Q , necessary to bring the theoretical prediction (||) into agreement 
with the given experimental value. The last line gives the present theoretical prediction. 



Reference 


Method 


r(o-Ps) [fis- 1 ] 


B 


Ann Arbor M 


Gas 


7.0514(14) 


338(36) 


Ann Arbor |J 


Vacuum 


7.0482(16) 


256(41) 


Tokyo 


SiC>2 powder 


7.0398(29) 


41(74) 


Tokyo (preliminary) || 


Si0 2 powder 


7.0398(15) 


41(37) 


QED prediction § 


7.0399 


44.87(26) 



The gas and vacuum results led to the suspicion that the two-loop QED correc- 
tions, which were not fully known, may be very large. The coefficient of (a/ V) 2 in 
the correction relative to the lowest-order decay rate is denoted by B Q . Its value, 
necessary to reconcile a given experimental result with the QED prediction, is given 
in the last column of Table IB. 

If such unusually enhanced effects existed, one would expect them to modify the 
p-Ps decay rate as well. Since p-Ps decay is much simpler than that of o-Ps, we 



2 



undertook a full two- loop QED study of this process jlClllTfl . 

The parapositronium decay rate into two photons agrees with predictions and 
had previously attracted less theoretical attention. However, it is also measured 
sufficiently precisely [ p"2|1 , 

r; x _ p Ps (gas) = 7990.9(1.7) /is" 1 , (1) 

to warrant a calculation of 0(a 2 s ) corrections. The prediction for this decay width 
can be parameterized as 
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Our aim was the evaluation of the second order non-logarithmic correction B p . It 
receives contributions from both soft and hard scales, B p = B s p oit + 5p ard . We found 



B s ° lt + 2vr 2 In - = — + 2tt 2 In 
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n hard 
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(3) 



and the final result is B p = 5.1(3). There is a significant cancelation between soft and 
hard pieces and the final result is almost eight times smaller than the magnitude of the 
finite constant in the hard scale contribution computed in dimensional regularization. 
Using the above result for B p we arrive at the following result for the decay rate: 



-ptheory 
1 p-Ps 



7989.64(2)/is~ 1 , 



which agrees very well with the measured value, Eq. (|T|). 

Most recently, the two-loop corrections were also computed for the o-Ps decay, 
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7.03994(l)//is, A = 10.286606(10). 



(4) 



The non-logarithmic two-loop coefficient is (with the small light-by-light contribution 
shown explicitly) B = 44.52(26) + 0.350(4) = 44.87(26) @,(Lj. For both o-Ps and 
p-Ps, the leading logs in the order ft 3 were found in |TJ]] and the next-to-leading 
logarithms were computed only recently |115l , l6"l ] l7]. 
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Can we claim that the positronium lifetime puzzle has been solved? It would 
certainly be very valuable to have another measurement of the o-Ps lifetime, especially 
that the powder measurement is somewhat controversial. An independent calculation 
of the theoretical prediction would also be useful, although it is unlikely that any 
unusually large effects are there to be uncovered. 



3 Expansion of bound-state energies in the constituent mass 
ratio 

A new approach to computing energy levels of a non-relativistic bound-state of 
two constituents with masses M and m, by a systematic expansion in powers of m/M, 
was described in [ P~8| . 

Simple atoms of experimental interest often consist of two particles (constituents), 
widely separated in mass. An extreme example is the hydrogen, where the ratio of 
the proton and electron masses is of the order of 2000. Smaller ratios characterize 



muonium, muonic hydrogen, and exotic hadronic atoms. The goal of |L8[ was to find 
a practical algorithm which would allow evaluation of the bound-state energy levels 
(in a given order of perturbation theory in a and Za) as an expansion in powers and 
logarithms of m/M with an arbitrary precision. 

That algorithm is useful in finding the so-called "hard-scale corrections". In the 
language of an effective field theory, this corresponds to determining the Wilson co- 
efficients of the short-distance operators, generated by virtual momenta much larger 
than the characteristic momenta of the atomic constituents. In practical terms, what 
is needed is the scattering amplitude of the two particles, with masses m and M, 
at the threshold, that is when the particles have vanishing velocities. The relevant 
Feynman diagrams depend on only the two mass scales m and M (since the external 
spatial momenta vanish). Since, however, at the two- loop level, which is of interest for 
the current theoretical studies, such integrals cannot in general be evaluated exactly, 
we need a method of expanding them in powers and logs of m/M. 

As is already well known in the theory of asymptotic expansions (for reviews and 



further references see e.g. ||19| , pO| , pT| ) , such expansion consists first of all in expanding 
the Feynman integrand in the small parameter m. In general, this gives rise to 
non-integrable singularities at small values of some momenta for which, before the 
expansion, the mass m provided a regulator. The difficulty in constructing a correct 
algorithm is to find the necessary counterterms and to evaluate the resulting new 
integrals. 

Such a procedure was carried out in |1| with the example of radiative-recoil 
corrections, such as the diagram shown in Fig. [I]. 

In the language of characteristic scales of the momenta in the counterterm inte- 
grals, there are 3 "regions". Fig. |2| depicts the Taylor expansion in m (@(a)) and the 
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Figure 1: An example of a forward-scattering radiative-recoil diagram. The bold line rep- 
resents the heavy constituent of the bound-state (e.g. proton if we consider hydrogen) and 
the thin line — the light one (an electron). 





(b) 




(c) 



(d) 



Figure 2: Elements of the expansion of the diagram in Fig. |l] in powers and logs of m/M 
(see text). Thick and thin solid lines denote propagators with M and m, respectively. Thin 
dashed lines are massless, and the thick dashed line denotes a static ("eikonal") propagator. 



3 counterterms. In the diagram ||(b), the momentum flowing through one of the light 
constituent lines is of the order of m, while the second loop momentum is O (M) (so 
that the propagators along its flow can be Taylor-expanded in m and in the corre- 
sponding small momentum). In 0(c) we have the same sizes of the loop momenta, but 
now the large momentum takes a different route to flow through the upper part of 
the diagram. The last contribution, 0(d), arises when both loop momenta are O (m). 
In this case, the heavy propagator l/(k 2 + 2k.P) can be expanded in k 2 , and becomes 
an "eikonal" propagator l/2k.P. 

This last contribution was least easy to evaluate and a procedure for computing 
the resulting integrals was described in ||18fl . It is interesting that the same integrals 
arise in very different problems, for example in certain corrections to the radiative 
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quark decays f23fl . 



4 Anomalous magnetic moment of a bound electron 



In [gj] the binding corrections to the gyromagnetic factor g e of an electron in 
hydrogen-like ions were studied. The interaction of an electron with an external 
magnetic field B is described by the potential V = —/j, ■ B. The electron magnetic 
moment [i is given by /x = g e j-s, with m, s = er/2, and g e denoting the mass, spin 
and the so-called gyromagnetic or Lande factor of the electron. 

If the electron is bound in a ground state of a hydrogen-like ion, g e becomes a 
function of the nuclear charge Z and its measurements |25,26,27j provide a sensitive 
test of the bound-state theory based on the QED. The theoretical prediction can be 
cast in the following form |28[ 



g e (Z) = g D + Ag vcc + Ag Illd . 



(5) 



The first term corresponds to the lowest order expansion in a/n and has been calcu- 
lated to all orders in Za |p9| , 



9b 
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(6) 



0((Za) 2 ^), where m N is the 



Ag rec denotes the recoil corrections |30| , Ag TCC 
nucleus mass. Further references to the studies of those effects can be found in p8 



Our main interest were the radiative corrections. They can be presented as an 
expansion in two parameters, Za and a/ir, 



Ag 



rad 



CP(Za) n+C^(Za) (°) +... 



(7) 



Powers of a/n correspond to electron-electron interactions, while Za governs binding 
effects due to electron interactions with the nucleus. The first coefficient function in 
(0); C^(Za), has been computed numerically to all orders in Za |TT],[52|. Its first 
two terms in the Za expansion are also known analytically |33|j34 



Cf\Z 



a 



l + -(Za) 2 + 0((Za) 4 
o v 



The main theoretical uncertainty for g e in light ions is, at present, connected with 
the unknown coefficient C in the next coefficient function, 



Cf\Za) = Cf)(0) [l + C'-(Za) 
C^(0) = -0.328 478 444 00... 



O [{Za) A 



(9) 
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At present, the most accurate experimental value of the bound electron gyromag- 
netic factor has been obtained f37jj38| with a hydrogen-like carbon ion 12 C 5+ (Z = 6), 

(10) 



g c {Z = 6; exp) = 2.001 041 596(5). 



The theoretical prediction is |39 



g c {Z = 6; theory) = 2.001 041 591(7) 



j a) 



where 70% of the error is caused by the unknown coefficient C of the (j^J {Zc 
effects in (^) (for carbon, higher powers of Za are assumed to be negligible). 

The purpose of our paper \24\ was to demonstrate that C — 1/6, in analogy to 
the corresponding coefficient in the lower order in a /it. In fact, we found that the 
coefficient of {Za) 2 is the same in all coefficient functions C^ 2n \Za), so that the 
theoretical prediction for Ag rad accurate up to {Za) 2 and exact in a/n reads 



Afi^ad — (<7f ree — 2) 



1 + 



{Zaf 
6 



:i2) 



where gi ree is the gyromagnetic factor of a free electron, presently known to O {{a/ir) 4 ) 
|f40|| (the same result had been obtained in a different way in |4T| ). With this result, 
the theoretical uncertainty in (|TTD is reduced from 7 ■ 10~ 9 to about 2 • 10 -9 . 



5 Conclusions 

Three examples of problems in the bound-state theory, solved since the last RAD- 
COR meeting, were summarized in this talk. Of course, other groups have made 
important progress in other aspects of bound-state physics, which was not reported 
here. I would like to mention two examples. Important new corrections to the hy- 
drogen Lamb shift were found in [[E^Q. In |42[| , an algorithm was constructed to 
compute a class of 3-loop propagator-type massive integrals (the so-called master 
integrals had been computed in [44)). 



Another class of recently studied problems is connected with the "velocity renor- 
malization group" f45j. This approach allows better understanding and at least par- 



tial resummation of large logarithms In - arising in the non-relativistic bound-state 
calculations. A review of the recent progress in this area can be found in . 



Hopefully, the recent significant progress in the theory will be followed by new 
measurements. It would be very important to re-measure the positronium hyper- 
fine splitting, to resolve the present (more than 3 standard deviations) discrepancy 
between the old measurements and the improved theory |P J4"?||i8| , |4*9| . Similarly im- 
portant would be an independent study of the ortho-positronium lifetime. Now that 
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the theory has in most cases taken control of the two-loop quantum effects, it is 
particularly exciting to test those predictions experimentally 
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